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1. Introduction and Summary 

The gauge/string correspondence and its numerous generalizations (see for a review) conjec- 
tures the equivalence between (super) string theory on a ten-dimensional space-time and some 
gauge theory living on the boundary of this space-time. The approach is based on taking the 
near horizon limit of the background produced by a large number of Dp-branes. On the field 
theory side the number of branes N c corresponds to the rank of the gauge group(s). There- 
fore the gauge/string correspondence provides a holographic description of the field theory 
in the planar limit. Large N c and supersymmetry contribute to the differences between the 
holographic models and the real world QCD. Yet, for many phenomenological applications the 
holographic approach can provide a valuable insight. In the holographic approach, in order to 
add a meson sector one has to include "flavor" D-branes ||, giving rise to the quark degrees 
of freedom. Such models, colloquially called holographic QCD, have been extensively studied 
over the last decade for various gravity backgrounds with different flavor D-branes. 

We choose the Klebanov-Strassler (KS) background as the starting point for our holo- 
graphic setup. This background is dual to a confining supersymmetric gauge theory which 
cascades to a Af = 1 SU(M) SYM in the deep IR. On the gravity side the KS solution is 
based on the 6d deformed conifold. The deformation parameter e is related to the dynamical 
scale A of the gauge theory. Apart from A, the real life QCD also includes additional massive 
parameters - the quark masses m q . In this paper we will focus on a holomorphic D7-brane 
previously studied in 0]. For the sake of simplicity, we will consider a single D7-brane. Thus 
we can use the quenched approximation ignoring the backreaction of the flavor brane on the 
gravity background. In general, the approximation is valid for N{ <^ N c . The holomorphicity 
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of the embedding implies that it breaks none of the background supersymmetries and the dual 
gauge theory with flavor is a Af = 1 SYM with fundamental matter. The embedding of the 
_D7-brane is parametrized by a single complex parameter p which is related to the quark mass. 
The preserved supersymmetry allows us to write down the precise form of the superpotential. 
It appears that apart from the quark mass term there is a term describing quartic interaction 
between the quarks and the KS bi-fundamentals. 

The low energy fields living on the D7-brane are the complex scalar that corresponds to the 
directions transversal to the brane and the 8d gauge field. From the Ad point of view the latter 
decomposes into a gauge field and another three real scalars. The gauge field is dual to the 
conserved U(l){ flavor symmetry current that gives rise to the vector mesons. All other fields 
are dual to scalar operators related to each other by supersymmetry. We review the meson 
multiplet structure later in the text. Although the calculation of the scalar meson spectrum 
currently appears to be a hard technical problem, the vector meson spectrum can be computed 
relatively easy for any value of p/e. We thus have a remarkable opportunity to investigate how 
the vector meson mass m depends on both the quark mass m q and the IR scale A, with the 
former fixed by p and the latter by e. 

It is worth recalling here that we are dealing with a supersymmetric version of QCD and 
as a result both A and m q (e and p) are complex parameters. Naively, for a fixed e, the meson 
masses are expected to grow as we increase |p|: the meson built of heavier quarks is supposed 
to be heavier. Surprisingly, we found numerically that this is not the case if the ratio x = p/e is 
real and smaller than a certain value, approximately equal to one. In this case only the lowest 
meson behaves as expected, i.e. its mass increases with x. At the same time the masses of all 
other modes decrease as x increases from to 1. Then, around x « 1 the masses reach their 
minimum and increase afterwards. This is an unusual behavior both because the meson mass 
is not a monotonic function of p and because it singles out the lowest mode. To confirm the 
result of numerics we use the WKB analysis reaching the same conclusion. 

In Nature the lightest vector mesons are the p-mesons built from the u and d quarks. These 
quarks are much lighter than the QCD scale. We compared the ratio of the lowest p-meson 
masses with the predictions of our model for p — and found a remarkably good agreement 
assuming that the p(2150) meson is a 3 3 Si mode. 1 

Next in the spectrum of vector mesons after the p's are the 0-mesons built from the heavier 
s and s quarks. These mesons, however, can mix the "pure" ss states with uu and dd. To 
compare with the holographic model we use putative pure "0" -mesons, whose masses can be 
derived from the spectrum of the K* and p-mesons. We found similarity between the non- 
monotonic behavior of the spectrum as a function of x in our model and the putative "0" -meson 
masses. Namely the first excited mode of the "0" -meson (that consists of the heavy s-quarks) 

lr The experimental status of the 3 3 Si mode is presently unclear. We discuss this issue later in the text. 
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is lighter than the first excited mode of the p-meson (that consists of the light u, <i-quarks). 
Remarkably, this applies only for the excited but not to the lowest mode, similarly to what we 
observed in our holographic model. 

These optimistic results should not be taken for granted since our model is only a distinct 
relative of the real-world QCD. Nevertheless, we believe that the observed similarities may not 
be a mere coincidence. It can rather suggest that the ratio of masses for certain light excitations 
in the QCD-like theories is quite robust and is not sensitive neither to the number of colors, nor 
to the explicit field content of the field theory. A similar phenomenon was previously observed 
for certain glueballs: the KS theory happens to predict with a reasonable accuracy the mass 
ratio for the lightest scalar glueballs of the pure SU (3) theory || . 

As another option, it may be the case that the inverse mass phenomenon is actually a result 
of some peculiarities of our model, not shared by QCD. Thus the quartic interaction between 
the quarks and the original fields of the KS model may provide a nice intuitive explanation for 
the different spectrum behavior for real and pure imaginary /x. This argument has also a clear 
geometrical counterpart. In terms of the dual geometry, the minimal distance between the tip 
of the deformed conifold and the lowest point on the D7-brane depends on x in an interesting 
way. When x has a non-zero phase the distance vanishes only for x — 0, while for real x the 
brane extends all the way to the tip as long as ^ x ^ 1. Hence we have different behavior for 
real and complex x. 

Organization of the paper 

In the next section we briefly review the Klebanov-Strassler (KS) model with the embedded 
_D7-brane. We also remind the reader the basics of the dual gauge theory. Next we discuss the 
_D7-brane world-volume fields and the meson multiplet structure. In Section we perform the 
calculation of the spectrum using the holographic approach. We compute the spectrum using 
the "shooting" technique and compare it with the results of the WKB approximation. Then 
we provide an interpretation of the spectrum behavior based on the peculiarities of our model. 
Section | is devoted to the comparison with the experimental data. 

2. The conifold based models: flavored and un-flavored 

The deformed conifold is a regular six dimensional non-compact CY manifold defined by the 
equation 

i=i 

The relevant 10c? supergravity solution is of the GKP type || with constant dilaton and a 
warped metric 

ds 2 w = h-^i^dx^dx 11 + h 1/2 (r)dsl . (2.2) 
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Here ds\ is the (deformed) conifold metric and the warp-factor h(r) depends on the transverse 
radial variable r defined by 



r 3 



En 2 - ( 2 - 3 ) 
i=i 

The minimal value of r is r = |e| 2//3 and in the deformed conifold case it is useful to introduce 
a new coordinate r defined by 

r 3 = |e| 2 cosh(r) , r > . (2.4) 

We refer the reader to the original paper [fj for the explicit form of the deformed conifold 
metric, fluxes and the warp factor. 

The field theory dual of the singular (e = 0) conifold background is a M = 1 Ad gauge 
theory at an IR conformal point, known as the Klebanov-Witten (KW) theory |7|. For the 
deformed conifold (e ^ 0) the gauge group is SU ((k + 1)M) x SU(kM) for integer k. The shift 
k k — 1 describes a single step in a cascade of Seiberg dualities, which reduces the gauge 
group down to SU (2M) x SU(M) in the deep IR. The R-symmetry is broken down to Z2 due 
to the formation of a gluino condensate (AA) ~ A 3 . On the gravity side the scale A is fixed by 
the conifold deformation parameter e. 

At the next step we consider a D7-brane embedded in the KS background. The D7-brane 
extends along the flat Minkowski space-time and a holomorphic non compact four cycle on the 
conifold E. The embedding we are interested in is 

Z4 = fi . (2.5) 

It was shown in [[| that the D7-brane wrapping E with trivial world- volume gauge field does 
not break the M = 1 SUSY of the original theory. Moreover, the 4-cycle wrapped by the brane 
is topologically trivial [§], §], thus the tadpole cancellation condition is satisfied and no anti 
-D7-brane is required in order cancel the net RR charge. On the gauge theory side it implies 
that the setup has no chiral symmetry. This is in contrast to the conifold based model of || Q 
(see also §|). 



The lowest value of the radial coordinate r m i n along the profile depends on the relative 
phase between \i and e. A simple computation shows that 

cosh (r min ) = |x 2 | + |l — x 2 \ , (2.6) 

where 

x = ^ . (2.7) 

This means that for real x we have r m i n = as long as \x\ ^ 1. On the other hand, if x has 
a non- zero imaginary part, then necessarily r m i n > 0. For fixed r the geometry of E is S 3 /Z 2 
which shrinks to S 2 at r = T m -, n . 
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The quiver diagram for the flavored Klebanov-Witten model is shown on Figure || and the 
superpotential is [[12], [L3|, [14[ 



W = \e ij e kl Tr {A^A^) - Tq{A x B x + A 2 B 2 - n)q - (qq) 2 , (2.8) 

where A\^, -81,2 are bi-fundamental matter fields. Thus the term 2fiqq is responsible for the 
quark mass. 

The right degrees of freedom for a confining theory are the di-baryon operators A and 
B constructed from the bi-fundamentals or mesonic operators of the form Tr (AiBj). One of 
the striking features of the KS model is that the superpotential retains its form under Seiberg 
duality. The KS solution describes an A -H- £>-symmetric point in the baryonic branch, which 
is evident, in particular, from the additional un-broken T, 2 symmetry of the background. 

Let us now address the map between various fields on the D7-brane and the gauge theory 
states. 

As we have already mentioned in the Introduction, the world-volume fields on the flavor 
_D7-brane are a complex scalar 5[i and a 8c? gauge field au- Throughout this paper we will 
adopt the a T = gauge. This leaves the Ad components and the legs along £ at finite 
t. The KS background describes a confining theory and only the colorless combinations of 
the quarks fields, like the mesons qq, appear in the spectrum. Both q and q are 4d chiral 
multiplets. The bottom component of the operator qq has dimension 3/2. This is evident 
from the superpotential ( |2.8| ). This operator corresponds to the combination a± + ia 2 . The top 
component J d 2 9 qq has dimension 5/2 and corresponds to the complex scalar describing the 
embedding. Furthermore, the world-volume gauge field with space-time indexes a M is dual to 
the £7(l)f vector current J M , and the remaining component of the world-volume gauge field 03 
is dual to the bottom component of the U(l){ current multiplet. 

In the holographic approach the 4g? mesons arise as the normalizable fluctuations of the D7- 
brane world-volume fields. The M = 1 SUSY arranges the particles into multiplets combining 
a massive vector with a real scalar or two real scalars together (we are not concerned with 



q Ai 




Figure 1: The quiver diagram of the flavored KW model. For the KS model the ranks of the two 
gauge groups are different and are not invariant under Seiberg duality. Yet, the structure of the quiver 
diagram survives the cascade of Seiberg dualities. 
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fermions here). There are no massless mesons in the setup in question. As was mentioned 
above the supersymmetry combines the vector mesons (the focus of this publication) with the 
scalars originating from the fluctuation of the world- volume gauge field 03. Although these 
particles have the same spectrum the equations they obey in the bulk are not the same. Rather 
they are related to each other through a supersymmetric quantum mechanics (SUSY QM) 
transformation [[TJJ. This situation is typical for the fluctuations that are dual to operators 
of different dimension. The equations for different particles in the multiplet coincide only if 
the corresponding dual operators are of the same dimension. In the present case the vector 
current has dimension 3 and the bottom component of the current multiplet has dimension 
2, and therefore only the related spectra will coincide but not the equations themselves. The 
SUSY QM transformation also relates the pair of scalars originating from the fluctuations of 
the world-volume gauge field ai j2 and the fluctuations of the embedding Re 5[i, Im S/i. Since 
a\ t 2 correspond to the operators of dimension 3/2 < 2, to find the spectrum one has to impose 
the "unusual" boundary conditions at infinity |TJ|. While the general asymptotic of is 
Cor -V2 _j_ Cir - 3 / 2 ; the proper boundary conditions for a\p are c\ — with arbitrary Cq- In 
practice this would mean that to obtain the spectrum from the original equations for a 12 in 
the bulk, using (for example) the shooting technique, may not be the optimal strategy. Rather 
one can study the equations for the superpartners Re S/jl, Im <5/i which yield the same spectrum. 
These fluctuations correspond to the operators of dimension 5/2 > 2 and although they have 
the same (up to an overall r~ 1//2 factor) asymptotic at infinity Cq + Cir -1 , they satisfy the 
conventional boundary conditions Co = with arbitrary c\. It happens that in the conformal 
KW case there is an accidental degeneracy between a\ and a 2 (as well as between Re dfM and 
Im 5fi) [|TIJ. It will be interesting to see if this degeneracy will be lifted in the KS case when 
the new parameter e is introduced. 



3. The spectrum calculation 

In this section we will calculate the spectrum of the vector mesons as a function of fi, e. These 
mesons are created by the U(l)t flavor current dual to the Minkowskian vector field living 
on the D7-brane. The five-dimensional action describing the lowest KK modes of this vector 
field has the form 

C ~ C(t)F^F^ + 2D(t)F„ t F? . (3.1) 

Here //, v are the Minkowski indexes and r is the radial direction along the conifold. We choose 
the gauge a T = and look for a solution of the form a M = v^e x ip(r) with v ■ k = and 

(D(t)iP'(t))' + m 2 C(r)^(r) = . (3.2) 

The four-dimensional meson mass m? = — k 2 emerges from the eigenvalue problem (|3.2|) . In 
general, to find the functions C{r) and D(t) could be a tiresome task. This is because both 
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the DBI and the CS parts of the D7-brane action contribute to the action (|3.1|) . The problem 
drastically simplifies for a SUSY D7-brane embedded into a constant dilaton, imaginary self- 
dual (ISD) background ||. In this case the invariant 2-form J 7 = ip*(B 2 ) + 2-Ka'dA, where ip* 
denotes the pull-back, satisfies the anti-self-duality condition 

T = — * 4 T , (3.3) 

which implies || 

C(T,/i,e)= VS-PW , and D(r, //, e) = g^y/g. (3.4) 

Here g stands for the un-warped induced metric on the 4-cycle £ wrapped by the D7-brane and 
*4 is the Hodge dual associated with g. The equation (|3.3|) follows directly from the K-symmetry 
condition, which in turn ensures that the embedding preserves the background supersymmetry. 
It was shown in [J3J that the -B-field in the original gauge of the KS solution fulfills the k- 
symmetry condition with zero world-volume gauge field A = 0. Hence in what follows we can 
substitute T with B assuming the pullback on £. 
As a result we find for our embedding 

g^Vg= |e| 4/3 (roA(r) ^ 2 ( ( r) ~ r)1/3 (sinh 2 (r)-2|x| 2 cosh(r)+x 2 + x 2 ) , (3.5) 



1 8/3 



^cosh(r) sinh(r) — r) 2 / 3 
8sinh 3 (r) 

2sinh 3 (Y) 



3(cosh(r) sinh(r) — r) 
2 1 / 3 |e|^_ 3 (rcosh(r) -sinh(r)) 2 , , _ _,, 

and the KS warp factor is 



x cosh(r) — x\ 2 + 

(sinh 2 (r) — 2|x| 2 cosh(r) + x 2 + x 2 ) 



PfB = 15 = Ixcosh(r) — x\ 

8 mi 2sinh 5 (r 



h(r) = ^—I(t), where J(r) = / dz *™"* 1 * -(sinh(2z) - 2zf/ z . (3.6) 

mi Jt sinh z 



~ 4/3 , , w s r° , zcoth^- 1 

'gb 

Here m g b is the scale mass of the gauge theory glueballs 



I I 2 / 3 

msb= 2v4^' (3 ' 7) 



and x is the complex dimensionless parameter defined above 



x = (3.8) 
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To derive the first expression in ( p.5|) one has to introduce a complete set of local coordinates on 
£, for example as it was done in The last two relations easily follow from the explicit form 
of the Kahler form and the NS flux on the deformed conifold as well as the relations between 
different forms on S flT8f . 

The mode i/j from ( |3.2|) is dual to an operator of dimension 3, a fact which is in agreement 
with its asymptotic behavior at infinity 2 

V>(r) « c + ci • e" 2r/3 . (3.9) 



The proper boundary conditions are Cq = 0, i.e. lim ip(r) = at infinity, and regularity at the 

origin r = r min , where ^ ~ 1 + g? 2 (t — r min ) 2 with a numerical coefficient d 2 = d 2 (x, m/m g b). 
These boundary conditions represent no challenge for numerical studies of the spectrum and the 
equation ( |3.2| ) can be solved with a standard shooting technique. What we actually calculate 
in this way is the ratio m/m g b between the meson and the glueball masses as a function of x. 
Since the glueball mass m gb is \i independent, changing x is equivalent to changing /i, while 
keeping e fixed. 

We solved ( |3.2| ) numerically for various values of x. Because the spectrum is invariant 
under x — > —x and x — > x* we focused only on the Re x, Im x > quadrant. In Table [I] we 
present our results 3 for x = and some real and purely imaginary values of x. As one can see, 
if x is real, for all modes, except for the lowest one, the meson masses decrease as x grows, 
until they reach their minimum value around x ~ 1. Then they start growing again. This type 
of behavior is not observed for the purely imaginary values of x, for which the masses grow 
monotonically with \x\. The results are presented graphically on Figures [| |3|. Our analysis 
seems then to imply that, for real x smaller than 1, the higher level mesons composed of heavier 
quarks are actually lighter than those composed of lighter quarks. 

This surprising behavior is confirmed by the WKB analysis, performed along the lines of 
[|l"9| (as reported in Appendix B of p0|). After putting equation (|3.2|) in the Schrodinger form 



<9 2 ( a/-D(t,/z, e)J e \ 

<9 2 0(r) - V{T,x,m/m gh )<f)(T) = , V{r,x,m/m gh ) = m 2 - ' 

y/D{T,n,e) D{T,fj,,e) 

(3.10) 

one can use the standard WKB formula 

fn — -jvr= / dt^J—V(t, x, m/m gh ) , for n ^ 1 . (3-11) 

\ / " T min 



2 Notice that at infinity r 3 sa \ |e| 2 e T 



3 The mass values for x — are different from those in [Q due to numerical errors therein. 
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X 






raj 


m gh 








n — 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 





1.487 


2.774 


4.077 


5.385 


6.696 


8.009 


0.15 


1.488 


2.773 


4.073 


5.379 


6.688 


7.999 


0.3 


1.490 


2.771 


4.061 


5.360 


6.661 


7.965 


0.45 


1.497 


2.767 


4.041 


5.326 


6.615 


7.907 


0.6 


1.504 


2.763 


4.014 


5.277 


6.547 


7.820 


0.75 


1.516 


2.763 


3.981 


5.211 


6.451 


7.697 


0.9 


1.531 


2.773 


3.957 


5.139 


6.326 


7.520 


1.2 


1.575 


2.880 


4.134 


5.395 


6.665 


7.941 


1.5 


1.629 


3.043 


4.427 


5.822 


7.226 


8.633 


1.8 


1.687 


3.208 


4.701 


6.203 


7.708 


9.216 


3 


1.909 


3.794 


5.626 


7.454 


9.279 


11.105 


7 


2.506 


5.249 


7.858 


10.450 


13.014 


15.584 


0.3 • i 


1.497 


2.802 


4.120 


5.443 


6.767 


8.094 


0.6 ■ % 


1.525 


2.878 


4.236 


5.599 


6.964 


8.328 


0.9 ■ % 


1.566 


2.985 


4.401 


5.820 


7.239 


8.659 


1.2 • i 


1.614 


3.108 


4.590 


6.073 


7.557 


9.041 


1.8 ■ i 


1.718 


3.373 


4.994 


6.613 


8.234 


9.851 


3 ■ i 


1.928 


3.886 


5.777 


7.660 


9.541 


11.421 


7 ■ i 


2.512 


5.279 


7.905 


10.503 


13.093 


15.679 



Table 1: The vector meson spectrum for real and imaginary x. The lowest value of the mass is 
indicated in bold. 



Here r* is the point where the potential vanishes, V(r*) = 0. For larger m this happens at 
larger r where V(r) can be approximated by 

Vir, x , mM „\ - (-£.) 2 (llzi) e -/3 + 0(o . (3 , 2 ) 

Thus r* is very large for large m and is x- independent. Hence we simply substitute it by infinity. 
In the large m limit we therefore have 

~ . y . n , where A(x) = / J 9(r, x)dr and 8(r, x) — -hr ^' . . (3.13) 
m gh A(ar) / Tmin m 2 gh D{r,fi,e) 

In Figure ^| we plot the meson masses calculated for real (left) and imaginary (right) x for the 
excited mode n = 7; the bottom solid line is the WKB result, while the upper dashed line is 
the result of the shooting. The strange behavior found for real x with the shooting technique 
is perfectly reproduced by the WKB approximation. 
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Figure 2: Meson masses m/m g b as functions of x = fi/e for n = 1 (left) and n = 7 (right). 




m/m gh 



n-1 \ xeiR 



ReOr) 




Figure 3: Meson masses for real (left) and imaginary (right) x for n = 7 from the WKB approximation 
(solid line) and the shooting technique calculation (dashed line). 



A similar behavior was found in the model studied in |TlJ (Figure 4). It is worth emphasizing 
that also in that model the bare quark mass can be tuned to be smaller than the dynamical 
generated scale. This is definitely not a common feature in holographyc models. For example, 
in the original Sakai-Sugimoto model |22| it is unclear how to tune the bare mass (the distance 
of the flavor branes from the tip of the cigar just reflecting variations of the constituent mass). 
In this sense, the KS background provides a privileged environment for the study of flavors 
with non-zero bare mass. 

We can derive the strange behavior of the excited (n > 1) modes for real x directly from the 
form of O(t) in ( |3.13| ). In the interval ^ x < 1 we have r m i n = and the function B(r, x) has 
the same asymptotic values at the origin 9(0, x) = 2- 5 / 3 3- 1 / 3 /(0) and infinity lim 6(r,x) = 0, 

T— >0O 

but it monotonically increases with x for any r in between zero and infinity as demonstrated 




Figure 4: The function 0(t, x) for x ^ 1 (left) and x > 1 (right). 



on the left graph on Figure ||. Hence the integral A(x) is a monotonically decreasing function 
with x, as expected. At x = 1 the asymptotic value of 0(r, x) at the origin r = jumps to 
0(0, 1) = 2~ 2 / 3 3~ 1//3 /(0) but A(x) is continuous. For x > 1 the asymptotic behavior of 0(t, x) 
changes drastically: it blows at r = r m ; n as (r — r m i n ) _1 . It is clear though that for x very 
close to 1 the behavior of 0(r, x) as a function of r will be the same as for x = 1 except for 
some very small vicinity of r = r min . Hence the integral A(x) is a continuous function of x 
as is expected from the physics of the problem. On the right graph on Figure [| we adjusted 
the minimal value of the coordinate r to the origin, f = r — r m ; n . It is clear now that as x 
increases the function 0(f,x) decreases at each f. Therefore A(x) is monotonically decreasing 
with x when x ~ 1. As a result we find that A(x) reaches its maximum, and the masses of the 
excited modes reach their minimum, at x = 1. This is close to the numerical results in Table [I], 
where the lowest masses of the excited modes appear around x = 0.9. It is worth mentioning 
here that the discussion in this paragraph is based on the WKB approximation which is usually 
not valid for small n, especially for n — 1, and so the monotonic behavior of the lowest mode 
cannot be addressed in this framework. 

The WKB analysis can also help us to see how the e dependence disappears in the KW 
limit when x ^> 1. Indeed in this case we expect the meson mass to depend only on /i and 
hence to be proportional to m g b|x| 2 / 3 . Such a behavior readily follows from ( |3.13| ). 6(r, x) goes 
as e~ 2r / 3 for large r and in order to find A(x) we need to integrate the square root of 0(r, x) 
from r m ; n to infinity. For large \x\ the parameter r m i n is also large and can be approximated 
by e Tmin 4|x| 2 as follows from (|2.6| ). Putting all together we find A(x) to be proportional to 
2 / 3 and therefore according to ( |3.13|) we get m n m m g b|x| 2//3 n as expected. 



\x 
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A possible explanation 

Before proceeding to the comparison with the known experimental data in the next section we 
would like to discuss the possibility that the strange behavior of the meson spectrum in our 
setup is a result of some intrinsic features of our model that are not present in QCD. 

Our proposal here is to attribute the observed phenomenon to the quartic interaction 
between the quarks and the bi-fundamental fields in the superpotential ( j2.8| ). Intuitively this 
interaction may modify the "effective" quark mass and hence the meson masses. Since we work 
in the quenched approximation, we might ignore the dynamics of the bi-fundamentals. Being 
more precise, in order to estimate the "effective" quark mass, we substitute the VEVs of Bi 
into ( |2.8|) . The VEVs are estimated as follows. 

The moduli space of the gauge theory is in one to one correspondence with the geometry of 
the deformed conifold. This fact is helpful in finding the superpotential (|2.8|) for the conformal 
KW theory. But the interpretation of fl2.8|) for the cascading KS theory is more tricky since 
the bi-fundamentals A4 and Bi are not the good degrees of freedom anymore. In the original 
paper || the "mesonic" fields = Tr(AiBj) were proposed as the low-energy fields in the 
dual gauge theory. It was later realized in J23j (see J24J] for the detailed discussion) that the 
theory is actually on its baryonic branch, My = and the low-energy fields are the di-baryon 
operators constructed from the bi-fundamental fields A,B~ \A\ ,\B\ . Although this step is 
not well-defined, we can evaluate the VEV of the original bi-fundamental fields A$ and Bi to 
be roughly of order y/e. This follows from the simple dimensional analysis. 

Thus, in order to estimate the "effective" quark mass, we should plug e instead of AiBj in 
(|2.8|) . The substitution gives roughly q(AiB t + A 2 B 2 — n)q — > q(e — fi)q, so that the effective 
quark mass is 

m e g ~ |e — fi\ — \e\ ■ |1 — x\ . (3-14) 

As we increase \i keeping e fixed, m e s strongly depends on the relative phase between /1 and 
e. In particular, for real x = /i/e the mass m e g decreases as we increase x all the way until 
we reach x — 1. This is closely related to the geometry of the D7-brane profile for real and 
imaginary fi as can be seen from the similarity between ( |3.14| ) and ( |2.6| ). Now, assuming that 
the usual dependence of meson mass on m e g (i.e. it increases when m e ff increases) we recover 
the unusual behavior of the excited modes for real < x < 1. An analogous explanation was 
suggested in ET] . 



In the light of this argument it is actually the behavior of the lowest mode in Table |I] that 
needs an explanation. Say, the lowest mode was precisely massless for /i = 0, i.e. it was a 
Goldstone boson of a certain spontaneously broken symmetry. Then for /i ^ the symmetry 
would be broken explicitly and the Goldstone boson would get a mass proportional to /i. Such a 
situation is realized in the model of || where the two massless modes correspond to the broken 
chiral and conformal symmetries. In our case there is no spontaneously broken symmetry at 
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fj, = 0. The behavior of the lowest mode, though, reminds of a pseudo Goldstone boson. This 
is similar to the situation with the pions in QCD, which are the pseudo Goldstone bosons of 
the chiral symmetry explicitly broken by the quark masses. 



4. Comparison to the experimental data 

In this section we compare our results with the spectrum of the p and mesons, i.e. the light 
quark 1 vector mesons in QCD. 4 Before proceeding let us recall that the correspondence 

e 2 / 3 ~A, // 2/3 ~m q , (4.1) 

is not rigorous. We just assume that the constant of proportionality is of the same order in the 
both cases. Similarly to e, the power of p here follows merely from dimensional analysis. 

Throughout this section we will use the standard n 2s+1 lj notations for the qq mesons | |25| | . 
Here n stands for the radial excitation number, s is either (anti-parallel quark spins) or 1 
(parallel quark spins), and I is the relative angular momentum of the quarks inside the meson. 
The latter is usually denoted by the conventional notations S, P, D, ... for I = 0,1,2,... 
Finally, \l — s\ ^ j ^ I + s is the total angular momentum of the state. An additional notation 
for the mesons is j PC , where P and C are the parity and the charge conjugation respectively. 
In terms of I and s one has P = (— and C = (— l) l+s , and so the S'-wave (I = 0) p-mesons, 
that we study in this paper, have 1 quantum numbers. 

1% errors for the p-meson masses 



The masses and widths of the first p mesons in units of MeV are J25 



p(770) 


775.49 ± 0.34 


(r = 


149.1 


p(1450) 


1465 ± 25 


(r = 


400), 


p(2150) 


2149 ± 17 


(T = 


359). 



(4.2) 



The first two are reported as 1 3 Si and 2 3 5*i modes; the latter is not firmly identified, but 
we report it for the reasons that will become apparent in a moment. We did not include the 
p(1700) because it is the 1 3 Di mode, and as such is not supposed to be captured by our analysis 
(the modes with higher spin in holographic models are parametrically heavier and correspond 
to stringy states; they are not captured by the supergravity approximation). We also ignored 
all heavier modes whose status is unclear. 



The heavy quark 1 mesons are not tightly bound in QCD, so the comparison with our values would be 
pointless. 
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We consider the two ratios 



^» ~ 1.88(9) , ~ 2.77(2) . (4.3) 

m p (770) ?™p(770) 



The p-mesons are mixtures of the u and d quarks. Hence the quark mass scale is m u ~ rrid 
(m u + md)/2 ~ 3.8 MeV. For the dynamical scale we can use A ~ m p ( 145 o) ~ 775 MeV as usual 
in the 7r-meson physics. 5 We will assume that 3.8/775 is the counterpart of the geometrical 
parameter x 2 / 3 . As we mentioned above the relation is not precise. Nevertheless, the basic point 
is that the resulting geometric ratio x is clearly very small, to the extent that using different 
prescriptions for x would not produce a sizable difference. 6 We are therefore allowed to use the 
x = results from Table [1] for the comparison 

^ „ 1.86(6) , ~ 2.74(2) . (4.4) 

m n= i m n=2 



The deviation of the values in ( |4.4| ) from the experimental data in (|4.3| ) are about 1% in both 
cases. In particular it is of the same order as the experimental error. It is now apparent why 
we included the p(2150) meson: the agreement of its mass with the mass of our third mode is 
so precise that we are led to believe that p(2150) is really the 3 3 Si mode. 

It is important to emphasize here that the status of the 3 3 Si mode is currently unclear 
(see for example p6|). Although an additional resonance(s) (with the same quantum numbers) 



between p(1450) and p(2150) has by now been observed by many groups, its nature and proper- 
ties, however, are far from being clear. This mode, sometimes called p(1900), with a predicted 
mass of about 1830 MeV, should exist if the spectrum of the light m q <C A mesons (like the 
ones we study) is consistent with the linear confinement, namely if m 2 ~ n. If indeed the linear 
relation m 2 ~ n is correct then the p(2150) meson will be the 4 3 S\ mode. This would reduce 
the precision of our results. 

The idea that m 2 grows as n can be supported by a simple semiclassical argument (see |27| 
for the basic review and the related AdS/CFT discussion). There are also some experimental 
indications that the masses of the light mesons satisfy the linear m 2 ~ n relations 



In the holographic models, however, m n usually grows like n 29 . In our case this is evident 



from the WKB approximation (|3.13| ). The possible exception being the lowest modes for which 
m 2 ~ n 2 does not necessarily hold. Hence the holographic models are not expected to describe 
the highly excited states with large n, but theoretically could give a reasonable prediction for 
the lowest states. 



5 A is supposed to set the scale for the meson resonances in QCD, modulo small corrections due to the finite 
quark masses, so a fair estimate is given by the first vector meson mass. 

6 For example, using A ~ 200 MeV we get x ~ 0.0026 but the resulting difference in the meson mass ratios is 
only of order 10 -3 . 
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Leaving aside the p(1900) issue, the surprising matching of the holographic calculation with 
the experimental meson masses is a priori unjustified. A list of reasons why our analysis may 
not be applicable includes the fact that the Klebanov-Strassler theory is supersymmetric, it 
contains many more light fields than QCD, has a large number of colors, and so on. 7 Given 
these fundamental differences, one possibility is that the remarkable agreement observed above 
happens by chance. Another explanation would be that, once the masses are normalized by the 
lightest one, these ratios are in fact robust and are not sensitive to the differences between the 
holographic model and QCD. Clearly, this can be true only for the few lightest modes before 
the different n dependence kicks in. Let us note here that a similar matching for the ratios of 
masses of the light(est) glueball states in the same holographic model and the pure glue SU(3) 
theory on the lattice was observed in f5|. 

Higher meson mass inversion 

The comparison with the mesons is more complicated, since the holographic calculation would 
correspond to the "pure" ss states, while the real mesons can be mixtures of the s with the 
u, d quarks. We will use a prescription from [[32|, [33|, [34| to give an estimate for the mass of a 



putative "pure" ss mesons (which we will denote as "0" -mesons in what follows 8 ). The mass of 
"0" is given by twice the mass of the corresponding K* meson (the light-strange vector meson) 
minus the mass of the p meson (the light-light) meson. 9 
The masses of the first two K* mesons are 

K* (892) : 896 ±0.25 (r = 50.3) , 
iT(1410) : 1414 ±15 [T = 232) , (4.5) 

and we obtain for the "0" 's 

"0(1020)" : 1015 ±0.5, 

"0(1680)" : 1363 ±33. (4.6) 

Note that the mass of the real 0(1020) is very close to the putative "0(1020)" one, consistently 
with the expectation of a very small mixing of the s~s state with the uu, dd for the ground state 
meson. Instead, the mass of the excited 0(1680) is quite larger than the one of the putative 
"0(1680)", a fact which agrees with a possible sizable mixing in this case. 



7 It is probably not even in the same universality class . The recent paper 1 31 discusses another source of 
discrepancy between the two theories, i.e. the effect on the p spectrum of the mixing of the vector and tensor 
quark bilinears in QCD, which is not captured by standard holographic models. 

8 As in [ p3[ the quotation marks are to remember that these are not the real particles. 

9 It is an implicit assumption here that the binding energy is more or less the same for both mesons. 
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^From these considerations we immediately obtain one of the main results of this note: 
the lowest mode of "0" in ( |4.6|) is heavier than that one of p in (|4.2| ) as expected, but the 
first excited mode of "0" is lighter than the one of the p-meson. This reproduces the behavior 
observed in Section [| for real x ^ 1. Indeed, one can see from Table 1 that the second (first 
excited) mode (n = 2) of our p meson, (i.e. essentially the x = entry equal to 2.774) is heavier 
than the second mode of our "0" mode (some entry for small x, e.g. for x — 0.15 the entry is 
equal to 2.773). 

A quick look at the measurement errors shows that they are not sufficient to account for 
this behavior in QCD. The difference in masses is not small, it is of order of 5 — 10% of the 
meson masses, and it is hard to attribute it to a mere coincidence. Of course, we are just 
considering putative "pure" ss mesons and, most importantly, the widths of such states are 
typically very large. Nevertheless, the effect described above is also quite large, so we believe 
that the qualitative behavior of our holographic model makes sense and it is actually not 
excluded that it could even be found in Nature. 

As an aside, let us add that if we want to compare the QCD "0" masses with the ones of the 
holographic model at hand we need to know the precise mapping of the geometric parameter 
p/e to the field theory value of m q /A. With all the reasonable choices, the discrepancy between 
our results and the QCD ones is larger than the one found for the p mesons and of the typical 
order for holographic computations applied toward phenomenology (5%-50%). Yet the crucial 
phenomena stressed above - the mass inversion of the higher meson modes i.e. that the higher 
mesons made of heavy quarks tend to be lighter than those made of lighter quarks, is not 
sensitive to the details of this map. 



Let us conclude this section by mentioning that in |20|, |35j the screening effects of light 
dynamical flavors on the spectrum of mesons composed of massless quarks in the Klebanov- 
Strassler theory has been analyzed. Along the same lines, it would be interesting to see the 
effects of light dynamical flavors on the results above and in particular on the observed behavior 
of the meson masses for real x < 1. 



Acknowledgements 

We are grateful to Francesco Bigazzi for collaboration on the early stages of this project. It is 
a pleasure to thank E. Klempt, J. Maldacena, D. Melnikov, A. Paredes and J. Soto for useful 
discussions, and J. Sonnenschein for his comments on the final version of the manuscript. We 
are grateful to E. Kiritsis for correspondence regarding the p-meson spectrum. 

A.D. thanks the theory group at Vrije University Brussel for hospitality while this work 
was initiated. The research of A.D. was supported by the Stanford Institute for Theoretical 
Physics, by the DOE grant DE-FG02-90ER40542, by the Monell Foundation, and in part by 
the grant RFBR 07-02-00878 and the Grant for Support of Scientific Schools NSh- 3035.2008.2. 



- 16 - 



The work of S. K. is supported by the European Commission Marie Curie Fellowship under 
the contract IEF-2008-237488. 

The research leading to the results in this paper has received funding from the Euro- 
pean Community's Seventh Framework Programme (FP7/2007-2013 under grant agreement 
n. 253534). This work is also supported by the FWO -Vlaanderen, project G. 0235. 05 and 
by the Federal Office for Scientific, Technical and Cultural Affairs through the Interuniversity 
Attraction Poles Programme (Belgian Science Policy) P6/11-P. 

A. L. C. would like to thank the Italian students, parents, teachers and scientists for their 
activity in support of public education and research. 



References 



F. Bigazzi, A. L. Cotrone, M. Petrini, and A. Zaffaroni, "Supergravity duals of supersymmetric 
four dimensional gauge theories," Riv. Nuovo Cim. 25N12 (2002) 1-70, fcep-th/0303191 . 



A. Karch and E. Katz, "Adding flavor to AdS/CFT," JEEP 06 (2002) 043, |hep-th/0205236 



I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: Duality 



cascades and chiSB-resolution of naked singularities," JEEP 08 (2000) 052, |hep-th/0007191 



S. Kuperstein, "Meson spectroscopy from holomorphic probes on the warped deformed 



conifold," JEEP 03 (2005) 014, hep-th/0411097 



A. Dymarsky, D. Melnikov, and A. Solovyov, "I-odd sector of the Klebanov-Strassler theory," 



JEEP 05 (2009) 105, |0810.5666 



S. B. Giddings, S. Kachru, and J. Polchinski, "Hierarchies from fluxes in string 
compactifications," Phys. Rev. D66 (2002) 106006, Eiep-th/0105097 . 



I. R. Klebanov and E. Witten, "Superconformal field theory on threebranes at a Calabi-Yau 



singularity," Nucl. Phys. B536 (1998) 199-218, [hep-th/9807080 



S. Kuperstein and J. Sonnenschein, "A New Holographic Model of Chiral Symmetry Breaking," 



JEEP 09 (2008) 012, 0807.2897 



A. Dymarsky, S. Kuperstein, and J. Sonnenschein, "Chiral Symmetry Breaking with non-SUSY 
D7-branes in ISD backgrounds," JEEP 08 (2009) 005, 



0904.0988. 



C. A. B. Bayona, H. Boschi-Filho, M. Ihl, and M. A. C. Torres, "Pion and Vector Meson Form 



Factors in the Kuperstein- Sonnenschein holographic model," JEEP 08 (2010) 122, |1006 . 2363 
M. Ihl, M. A. C. Torres, H. Boschi-Filho, and C. A. B. Bayona, "Scalar and vector meson 



spectra of flavor chiral symmetry breaking in the Klebanov-Strassler background," 1010.0993. 



-17- 



[12] 

[13] 

[14] 

[15] 
[16] 

[17] 
[18] 

[19] 

[20] 

[21] 
[22] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 



P. Ouyang, "Holomorphic D7-branes and flavored N = 1 gauge theories," Nucl. Phys. B699 



(2004) 207-225, |hep-th/0311084 . 

F. Bigazzi, A. L. Cotrone, and A. Paredes, "Klebanov-Witten theory with massive dynamical 



flavors," JHEP 09 (2008) 048, 0807.0298 



F. Bigazzi, A. L. Cotrone, A. Paredes, and A. Ramallo, "Non chiral dynamical flavors and 



screening on the conifold," Fortsch. Phys. 57 (2009) 514-520, |0810.5220 . 
F. Benini and A. Dymarsky; in progress. 

I. R. Klebanov and E. Witten, "AdS/CFT correspondence and symmetry breaking," Nucl. 



Phys. B556 (1999) 89-114, |hep-th/9905104 . 

F. Benini et. al, "Holographic Gauge Mediation," JHEP 12 (2009) 031, |0903.0619| . 

A. Dymarsky, "Flavor brane on the baryonic branch of moduli space," JHEP 03 (2010) 067, 



0909.3083 



J. G. Russo and K. Sfetsos, "Rotating D3 branes and QCD in three dimensions," Adv. Theor. 
Math. Phys. 3 (1999) 131-146, frep-th/990"i~056 . 



F. Bigazzi, A. L. Cotrone, A. Paredes, and A. V. Ramallo, "Screening effects on meson masses 



from holography," JHEP 05 (2009) 034, |0903.4747 . 

A. Paredes, "On unquenched N = 2 holographic flavor," JHEP 12 (2006) 032, [hep-th/061027C| . 
T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. Theor. 



Phys. 113 (2005) 843-882, |hep-th/0412141 



O. Aharony, "A note on the holographic interpretation of string theory backgrounds with 
varying flux," JHEP 03 (2001) 012, fcep-th/OloToil . 



A. Dymarsky, I. R. Klebanov, and N. Seiberg, "On the moduli space of the cascading SU(M+p) 



x SU(p) gauge theory," JHEP 01 (2006) 155, hep-th/0511254 



Particle Data Group Collaboration, S. Eidelman et. al, "Review of particle physics," Phys. 
Lett. B592 (2004) 1. 

E. Klempt and A. Zaitsev, "Glueballs, Hybrids, Multiquarks. Experimental facts versus QCD 
inspired concepts," Phys. Rept. 454 (2007) 1-202, |0708.4016| . 



A. Karch, E. Katz, D. T. Son, and M. A. Stephanov, "Linear Confinement and AdS/QCD," 
Phys. Rev. D74 (2006) 015005, |hep-ph/0602229| . 



A. V. Anisovich, V. V. Anisovich, and A. V. Sarantsev, "Systematics of q anti-q states in the 
(n,M**2) and (J,M**2) planes," Phys. Rev. D62 (2000) 051502, lhep-ph/00031ll . 



- 18 - 



[29] M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, "Meson spectroscopy in AdS/CFT 



with flavour," JEEP 07 (2003) 049, |hep-th/0304032 



[30] S. S. Gubser, C. P. Herzog, and I. R. Klebanov, "Symmetry breaking and axionic strings in the 
warped deformed conifold," JEEP 09 (2004) 036, [hep-th/0405282 . 



[31] S. K. Domokos, J. A. Harvey, and A. B. Royston, "Completing the framework of AdS/QCD: 



h\jb\ mesons and excited omega/rho's," 1101.3315 



[32] C. R. Allton, V. Gimenez, L. Giusti, and F. Rapuano, "Light quenched hadron spectrum and 



decay constants on different lattices," Nucl. Phys. B489 (1997) 427-452, [hep-l.at/9611021 



[33] I. Iatrakis, E. Kiritsis, and A. Paredes, "An AdS/QCD model from Sen's tachyon action," Phys. 



Rev. D81 (2010) 115004, |1003.2377 



[34] I. Iatrakis, E. Kiritsis, and A. Paredes, "An AdS/QCD model from tachyon condensation: II," 
JEEP 11 (2010) 123, 



1010.1364 



[35] F. Bigazzi, A. L. Cotrone, A. Paredes, and A. V. Ramallo, "The Klebanov-Strassler model with 
massive dynamical flavors," JEEP 03 (2009) 153, |0812.3399. 



- 19 - 



